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ABSTRACT 

The atmospheres of close-in extrasolar planets absorb most of the incident stellar radiation, advect 
this energy, then reradiate photons in preferential directions. Those photons carry away momentum, 
applying a force on the planet. Here we evaluate the resulting secular changes to the orbit, known 
as the Yarkovsky effect. For known transiting planets, typical fractional changes in semi-major axis 
are about 1% over their lifetime, but could be up to ~ 5% for close-in planets like OGLE-TR-56b 
or inflated planets like TrES-4. We discuss the origin of the correlation between semi-major axis 
and surface gravity of transiting planets in terms of various physical processes, finding that radiative 
thrusters are too weak by about a factor of 10 to establish the lower boundary that causes the 
correlation. 

Subject headings: celestial mechanics — planetary systems — planets, individual (HD209458b, 
HD189733b, HAT-P-2b, OGLE-TR-56b, TrES-4) 
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1. INTRODUCTION 

The radiative environment of a close-in extrasolar 
planet is intense. At the canonical hot- Jupiter dis- 
tance of 0.05 AU around a solar-type star, the incoming 
flux is a factor of 10 4 larger than that experienced by 
Jupiter itself. The starlight incident on a hot Jupiter is 
mostly absorbed, as very strong upper limits on reflected 
starlight attest: the geometric albedo of HD209458b 
at visi ble wavelengths is less than 0.17 at 99.5% confi- 
dence (jRowe et alJ l2007h . much smaller than Jupiter's 
value of 0.5. As the bulk of the planet is expected 
to be tidally locked, the time-averaged radiation field 
is strongly anisotropic. This heat input drives winds 
around the planet that are of planetary scale. 

The atmospheric dynamics on hot Jupiters is a hot 
topic, for both theory and observation. A consistent the- 
oretical prediction of three-dimensional models is that, in 
a steady state, the hottest portion of the planet does not 
lie at the substellar point, but is shifted east ward by a 
superrotating wind (Showman fc Guillotjr2 002). The the- 
oretical interpretatio n of this feature of the simulations 
is still uncertain, but lShowman et ail (|2007l ) suggest the 
day-night temperature contrast triggers eddies that drive 
eastward angular momentum to the equator. Observa- 
tions are starting to reveal the existence and extent of the 
phase shift. A hotspot displaced to the east by ~ 30° was 
recently observed at 8 Jim on the planet HP 18 9733b by 
iKnutson et "all (|2007[ ). Harrington et all (|2006l ) observed 
the phase curve of v And b to have a best-fit westward 
shift of ~ 11°, but the data are roughly consistent with 
no shift. 

The orbital mechanics of hot Jupiters depend on the 
radiation from their atmospheres. Because light carries 
momentum, it exerts a force as it is absorbed or emitted 
anisotropically. Absorbing the incident light pushes the 
planet away from the star. Since the amount of radiation 
intercepted is proportional to the inverse square of sep- 
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aration, the effect is the same as a smaller stellar mass, 
but the orbit remains Keplerian. On the other hand, if 
the radiation from the planet is preferentially tangential 
to the orbit, its semi-major axis can increase or decrease. 
If it is true that the hotspot moves eastward on average, 
the reradiated light acts as a thruster back along the path 
of the planet, adding angular momentum to the planet, 
which promotes it to a higher orbit. 

Thermal radiation acting as a thruster which changes 
the orbit is well-known in asteroid dynamics as the 
Yarkovsky effect (|OpikHl95lHBottke et aDl2006h . In that 
case, sunlight heats up the day side, but the asteroid's 
spin is not synchronous with its orbit, and rotation dis- 
places the energy before it is reemitted. This effect causes 
substantial orbital changes over the age of the solar sys- 
tem only for small (< 10 km) asteroids due to their large 
cross-section per unit mass. Relative to a 10 km ra- 
dius asteroid in the main belt of the solar system, hot 
Jupiters are ~ 10 12 times more massive. However, they 
are also ~ 10 2 times closer to the star and have a ra- 
dius ~ 10 4 times larger; i.e., they intercept 10 12 times 
the luminosity. Therefore, if the reradiation occurs with 
similar phase angle and contrast, the fractional orbital 
changes of planets close to other stars will be the same 
as that of asteroids of the solar system. 

The plan of this Letter is as follows. In $2]we show the 
rate of change of orbital parameters according to celestial 
mechanics. In fj3]we calculate these effects for two model 
planetary atmospheres and an observed planet. In i|4]we 
compare the observed distribution of transiting planets 
to limits due to various physical processes that tend to 
move or destroy them. Finally, S}5] closes with a brief 
discussion and call for further study. 

2. EFFECT ON THE ORBIT 

Here we compile for mulae that give the effects of a force 
on an orbit, following lMurrav fc PermottJ (|1999( ), section 
2.9. The force imparted by the radiation can be decom- 
posed into radial, transverse, and normal components: 
F = Rr + TO + Nz. 
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The rate of change of semi-major axis is: 
2a 3 / 2 



m p ^G(AU + m p )(l - e 2 ) 



[Re sin/ + T(l + ecos/)], 
(!) 

where / is the true anomaly, e is eccentricity, and M+ 
and uip are the stellar and planetary masses. The thrust 
is produced by the dipole of the reradiated flux distribu- 
tion, which we characterize by its fractional amplitude 
fd, so that the flux distribution has an angular depen- 
dence cx (1 + fd cos ip), where ip is measured from a pole 
through the peak of the dipole. The incoming starlight 
is Doppler-boosted in the frame of the planet, and upon 
absorption causes angular momentum loss (Poynting- 
Robertson drag), but that effect is v/c ~ 5 x 10~ 4 weaker, 
where v is the orbital velocity and c is the speed of light. 
The force due to the dipole component is 



8c V r / 



(2) 



where £+ is the star's luminosity, R p is the planet's ra- 
dius as measured by photometry during transit, r = 
a(l — e 2 )(l + ecos/) _1 is the distance between the planet 
and star, and A is the planet's albedo. (We assume that 
the albedo does not vary over the surface of the planet; 
if it does, there is additional opportunity for momentum 
transfer between starlight and the orbit.) We assume 
atmospheric dynamics displace the dipole a longitudinal 
angle Ad east of the substellar point, so the components 
of this force are: R = F cos Ad, T = F sin Ad, and N = 0. 
For an order-of-magnitude estimate, let us assume Ad and 
fd are constant. There may be an anti-correlation in the 
magnitudes of Ad and fd as high-altitude opacity sources 
may cause planets with high effective t emperature to 
reradiate the energy before it is advected (jFortnev et al.l 
2007a|). Averaging equation (P) over time, we have: 



A) /d sin Adi? 2 



4cG 1 /2(il/^ + m p ) 1 /2 mpa i/2(i _ e 2) ■ 
Therefore the timescale for semi-major axis change is: 
r a = a I a 



(3) 



=5.5 x io" J^M-V^^))- 1 /!^) 

■ y V 0.5 J V C & J \mjR- 2 ) 
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a \3/ 2 /M* + m p \V 



0.05AU 
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where reference values are for Jupiter and the Sun. 
Therefore the timescale of orbital change for a typical 
hot Jupiter is about 100 times its age. 
The rate of change in eccentrity is: 



e = — J " ( ! P J flsin/ + T(cos/ + cosg)], (5) 
m p y G(M* + m p ) 

where £ is the eccentric anomaly. Now, using the same 
force and assumptions as above, we have: 



r e = (dlne/dt) 1 = r a 



2e 2 



l_ e 2_( 1 _ e 2)3/2 



4r Q , (6) 



where the final approximation is first-order in ec- 
centricity. The Yarkovsky effect causes eccentricity 
to grow if the hotspot is shifted to the east (e.g., 



iSpitale fc Gr"eenbe ra 2002), on a timescale similar to that 
of semi-major axis growth. The assumptions of constant 
fd and Ad will break down at large eccentricity, which 
will be shown in fj3j However, given that the Yarkovsky 
effect will only be large for small values of R p /a, tidal 
dissipation in the planet should be sufficient to damp 
the eccentricity faster than the Yarkovsky effect excites 
it. 

The rate of change in orbital inclination is: / oc 
N. Above we assumed that the reradiation field is 
anisotropic, but is still symmetric above and below the 
plane of the orbit. In reality, the normal component of 
the force, A, fluctuates with fluctuating eddies. How- 
ever, this is a random walk with extremely small step 
size, so it is unlikely to get very far, even if the normal- 
ization of the rate were large. 

3. CALCULATION OF THE FORCE 

In this section the rates of orbital evolution are com- 
puted for three planets, taking the flux distribution for 
HD 209458b and HAT-P-2b from atmospheric models 
and for HD 189733b from infrared observations. 

Atmospheric models readily give the effective temper- 
ature distribution T c ff of the planet, each patch of which 
radiates into a hemisphere. The thrust from a patch is 
SF = —f'crT^fSA/c. Here 5A is the outward normal 
vector to a patch of area 5 A, a is the Stefan-Boltzmann 
constant, and the factor /' accounts for the fact that the 
radiation is not all emitted normal to the surface: if it 
were, /' would be 1. If treated as a true blackbody, the 
radiation from a patch uniformly illuminates the hemi- 
sphere, and f' = h- On the other hand, a limb-darkened 
atmosphere preferentially emits straight up rather than 
sideways, so /' is between ^ and 1. For instance, in 
a plane-parallel gray atmosphere with intensity propor- 
tional to (fi + 2/3), where u is th e cosine of the angle 
to the zenith (e.g., Mihalasl ll978l ). this factor becomes 
/' = 4/7. Integrating SF over the surface of the planet 
gives the net force. Finally, integrating equations ([T]) and 
([5]) over an orbit gives the secular effect of this force on 
the orbital elements. 

For the planet HD209458b, ICooper fc Showman! 
(2006) carr ied out a detailed a tmospheric circulation 
model and iFortnev et al.l ((2006) ran radiative transfer 
calculations to predic t phase curves. For th e equilibrium- 
chemistry case from IFortnev et al.l (|2006h , we compute 
the dipole component of the flux from the planet (based 
on the hemisphere-averaged T e ff as a function of orbital 
phase) to have an amplitude of 4.9 x 10 28 erg s , with its 
maximum shifted an angle 29° eastward. The timescale 
for semimajor axis change is then 3.8 x 10 11 yr; for a 
fractional change of Aa/a = 0.3% over t he age of the 
syste m. (All age estimates are taken from iTorres et all 
120081 ') The non-equilibrium chemistry case has a darker 
night side, so the dipole component is bigger, enhancing 
t he Yarkovsky effect by abou t 20%. 

lLangton fc Laughlinl (|2007f ) have computed theoreti- 
cal atmosphere models for a number of eccentric hot 
Jupiters. Here we analyze their computation for the 
eccentric planet HAT-P-2b, choosing it because it (1) 
has a rather large eccentricity (0.5) which complements 
the analysis in <J3] and (2) is orbiting a bright star 
and will therefore receive detailed observations. Sec 
lLangton fc Laughlinl (|2007l ) for the details of the model; 
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■ HD147506b = HAT-P-2b 
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Fig. 1. — The force vectors due to radiation asymmetry as a 
function of orbital position for the planet HAT-P-2b, calculate d 
from the atmospheric simulations of Langton & Laughlin (2007). 
Only the vector components in the plane of the orbit are shown. 



we additionally assume is that there is no limb darkening. 
The model was run for 19.5 orbits to reach dynamical 
equilibrium, after which we computed the force as a func- 
tion of time over a single orbit. Figure Q] shows that force 
as a function of orbital position; neither A<2 nor fd are 
constant. This force causes growth in both semi-major 
axis and eccentricity with timescales of 7.6 x 10 12 yr and 
1.1 x 10 14 yr, respectively, considerably longer than the 
estimated age of 2.6 x 10 9 yr. The rather large eccentric- 
ity causes the force to lose coherence over the orbit, and 

T e > T a . 

The planet HD189733b has been observed for about 
60% of its orbit, and the 8 /xm-brightncss of the planet 
as a function of longitude has been evaluated in 12 long i- 
tudinal bins of constant brightness (|Knutson et alj [2007) . 
This 8 [Mm "map" needs to be converted to a bolometric 
flux "map" to compute the Yarkovsky effect. Most of the 
flux comes out within a factor of two in wavelength from 
8 /zm, so we assume these two maps would be roughly 
the same, although scaled by the bolometric flux ratio of 
the planet and the star, rather than that ratio at 8 /xm. 
We assume an albedo of 0.1 and no limb darkening. The 
net force due to reradiation is 1.7 x 10 19 g cm s -2 in 
a direction mostly outward from the star, but tipped in 
the planet's orbital direction by A^ = 6°. This shift is 
smaller than the eastward displacement of the hottest re- 
gion of the planet (~ 30°) because the coldest region is 
in the same hemisphere, lessening the component of the 
force along the orbit. The timescale for semi-major axis 
change is 3.0 x 10 12 yr, so over 5 Gyr, this thrust pro- 
duces a 0.17 % increase in the semi-major axis. Stellar 
models suggest the star is at least a few Gyr old. 

In the future, both observational and theoretical work 
on very close planets should compute and report the mag- 
nitude of the radiative thruster effect by the methods 
illustrated here. 




0.01 0.10 
semi — major axis [AU] 

Fig. 2. — The surface gravity versus semi-major axis for known 
transiting planets. Various theoretical mechanisms disfavor the 
existence of planets on the hashed sid es of each label l ed line. The 
dashed lines are gas-giant models by Fortncy ct al. (20073); see 
text. 

4. GRAVITY/SEMI-MAJOR AXIS CORRELATION 

The scaling of equation (g]) suggests that the rel- 
evant property o f the planet is its surface gravity. 
iSouthworth et ail (|2007f ) have shown that an anti- 
correlation exists between surface gravit y and semi- 
major axis for transitin g planets (see also iMazeh et all 
120051: iTorres et all 120081) . Figure is a plot of the rela- 
tion updated with newly discovered planets. The sample 
has a sloping lower boundary, and apart from three out- 
liers, a clear trend. Lines mark where typical timescales 
of theoretical effects equal 5 Gyr, with a hashed area to 
the "forbidden" side of each constraint. 

The Yarkovsky effect, according to equation (|4]) using 
the nominal values of all parameters, is represented by 
t he line labeled "rad iative thrusters" . 

iRasio et al.l (|1996f) raised concerns that very close plan- 
ets would be pulled in towards the star due to tidal dis- 
sipation in the star's enve lope. The timescale for tidal 
decay is (jRasio et al.| [l996) : 

where ft has a very uncertain value depending on the 
tidal dissipation efficiency. This constraint for a R p = R,j 
planet orbiting a solar-type star is plotted in Figure 
Although the radiative thruster effect with an eastward- 
shifted hotspot tends to push the planet out, tides have 
a much stronger radial dependence, therefore the equi- 
librium between the two effects is unstable. 

The tidal gravity of the star limits how close a planet 
of a given radius can orbit without being disrupted; this 
"Roche limit" is shown in Figu re [2] for a lRj planet on 
a circular orbit. iFord fc Rasiol (|2006l ) have pointed out 
that the hot Jupiters obey an inner boundary at about 
twice their Roche limit, which implicates tidal circular- 
ization from large eccentricity and semi-major axis. If 
that mechanism delivers most hot Jupiters to their cur- 
rent orbits, then it also strongly affects their structure, 
especially for the least massive ones. Below a few J upiter 
masses, a tidal inflation instability (|Gu et al.ll2003f ) can 
result, which can strip the gaseous envelope and push out 
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planets to the Roche limit computed with the severely 

inflated radius. 

Models of atmospheric escape are very uncertain, but 

they suggest a hydrodynamic wind is being driven out of 

the planet's R oche lobe by extreme ultraviolet photons 

from the star. lYelld (|2004h found by use of a chemical 

reaction network that matter can be lost at a rate of 

4.7 x 10 10 g s" 1 for HD209458b. We may scale this rate 

to planets of different mass and radius, although several 

assumptions are necessary. First, we assume that the 

optically-thick area of the absorbing species in the exo- 

sphere is A oc R p [probably A w (3i? p ) 2 ]. Assuming that 

the escaping wind uses its energy to climb out of the 

planet's potential well, rn p oc ARpTn~ . Combining with 

the results for HD209458b gives the following timescale 

for mass loss: 

m p 
r m = ~~~ 

TO p 

Here the power on a is given by Figure 7 of lYdia (j200l . 
and is close to 2, the scaling due to a fixed cross section 
to XUV photons. This constraint is plotted in the figure 
as well, for R p — Rj. We note that the stellar XUV 
luminosity is strongly dependent on the spectral type of 
the star, which should be taken into account for evap- 
oration (jLecavelier Pes Etangd [2007) . The "?" symbol 
on Figure [5] means that the shape and normalization of 
the plotted limit is still widely debated. 

Also plotted in Figure [2] are g as giant models of thre e 
different masses from Table 4 of IFortnev et all (|2007b[ K 
The models are of solar metallicity gas with no rocky 
core, and have been contracting for 4.5 Gyr, while ir- 
radiated by a solar-type star (which produces the semi- 
major axis dependence). If close-in terrestrial planets 
can hold on to a thick atmosphere they may have sim- 
ilar atmospheric dynamics as gas giants, along with the 
Yarkovsky effect. The models of IFortnev et all (|2007bl ) 
for Earth-like compositions and masses (0.1, 1.0, 10)Af® 
have surface gravities of (0.16, 0.39, 1.22) in Jupiter units, 
so the orbits of terrestrial planets soon detectable by Ke- 
pler probably have not moved much by this effect. Such 
planets without an atmosphere will probably also avoid 
the asteroid-like Yarkovsky effect (due to spin plus ther- 



mal inertia) because of tidal spin synchronization. 

5. DISCUSSION 

The goal of this paper is to alert the exoplanet com- 
munity to a dynamical effect that is very important for 
understanding the orbits of small asteroids in the solar 
system and to suggest it may modify the orbits of the 
closest-in extrasolar planets. 

The results of $3] give long timescales for semi-major 
axis and eccentricity change for known hot Jupiters that 
have been studied extensively, much longer than their 
estimated ages, or indeed, the age of the universe. Of 
the observed planets, the biggest fractional semi-major 
axis change is for OGLE-TR-56 and TrES-4, which may 
have migrated by ~ 5 % during their lifetimes, ac- 
cor ding to equa ti on [^] using the parameters determined 
by iTorres et alj ([2008). The scalings of the radiative 
thruster effect imply that it can make a substantial 
change for the most extreme extrasolar planets, so at- 
mospheric circulation studies should be undertaken for 
such planets. 

Even a few percent orbital change could enhance the 
chance of success of current attempts to discover ad- 
ditional planets in mean-mot i on re s onance by transit- 
timin g variations <|Agol et al.1 120051: iHolman fc Murray! 
2005). The orbits of most known transiting planets 
have likely been modified by tidal dissipation in the 
planet, which tends to disrup t mean-motion resonance 
(|Terquem fc Papaloizoul feOOTf h However, the radiative 
thruster effect could oppose this motion, keeping planets 
in resonance, where transit timing can probe to much 
smaller masses. 

Finally, we emphasize that the remaining theoretical 
uncertainties will soon be settled by observations of plan- 
etary atmospheres, which will determine the importance 
of the Yarkovsky effect in the context of extrasolar plan- 
ets. 
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cussions, and A. Showman, S. Tremaine, M. Holman, and 
D. Sasselov for comments. These ideas were developed 
during a Michelson Fellowship supported by the National 
Aeronautics and Space Administration and administered 
by the Michelson Science Center. 
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